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We consider the dynamics of a thin hquid film in the attractive substrate potential and under 
the action of a longitudinal or a tilted vibration. Using a multiscale technique we split the film 
motion into the oscillatory and the averaged parts. The frequency of the vibration is assumed 
high enough for the inertial effects to become essential for the oscillatory motion. Applying the 
lubrication approximation for the averaged motion we obtain the amplitude equation, which includes 
contributions from gravity, van der Waals attraction, surface tension, and the vibration. We show 
that the longitudinal vibration leads to destabilization of the initially planar film. Stable solutions 
corresponding to the deflected free surface are possible in this case. Linear analysis in the case of 
tilted vibration shows that either stabilization or destabilization are possible. Stabilization of the 
dewetting film by mechanical action (i.e., the vibration) was first reported by us in PRE 77, 036320 
(2008). This effect may be important for applications. 

PACS numbers; 47.15.gm, 47.20.Ma, 68.08.Bc 



I. INTRODUCTION 

Dynamics of thin liquid films was extensively studied 
during the last decade both experimentally and theoret- 
ically. The importance of such studies is emphasized by 
the needs of modern nano and microfluidic technologies, 
which commonly employ films in the 100 — lOOOA thick- 
ness range. Reviews focusing on different subfields of 
research include Refs. [1,, ^ ^3,], as well as twelve reviews 
focusing on wetting in the recent volume [Annu. Rev. 
Mater. Res. 38, (2008)]. 

It is well-known that very thin liquid films tend to 
dewet from the substrate (rupture). The primary cause 
for dewetting is the attractive van der Waals interaction 
of the film and the substrate. Loss of stability and rup- 
ture of liquid sheets is often undesirable and may lead to 
a technological or manufacturing process failure. Thus 
understanding dewetting and finding means to control it 
are the important and challenging problems. 

One of the frequently used methods for controlling the 
fluid flow on small-to-large spatial scales is the applica- 
tion of the high frequency vibration [1, @. Several 
phenomena may emerge when such vibration is applied, 
such as the oscillatory (pulsatile) fluid motion (Faraday 
instability) and the time-averaged fluid motion. Analy- 
ses of the pulsatile motion of the liquid layer and thin 
drops are carried out in Ref s. M , 0, for the transversal 
vibration and in Refs. [gI. [qI. [ToI| for the longitudinal one. 

The standard high-frequency approximation is based 
upon the assumption of the vibration frequency so large, 
that the viscosity is important only in a thin boundary 
layer near the rigid wall 0, S • This approximation works 
well for macroscopic films [6, [ll|, [l^], but it fails for the 
thin films. Nonetheless, we have recently demonstrated 
that a hierarchy of typical times allows for the aver- 
aged description in a thin film system. Instead of using 



the standard high-frequency approximation, we assume 
that the vibration period 2tt/uj is (i) of the order of the 
characteristic time of the transversal transfer of the mo- 
mentum, Hq/i/, and (ii) is small compared to the typi- 
cal "horizontal time", L^/i/. (Here Hq is the mean film 
thickness, is the kinematic viscosity, and L is the typi- 
cal horizontal scale. L H for the thin film.) 

Reference [11] develops the averaged description for 
the case of the vertical vibration of the substrate. We 
show that the influence of the vibration is flnite if the 
amplitude is large in comparison with Hq- In this case the 
vibration is the efficient way to stabilize the film against 
the van der Waals rupture. In this paper the approach 
of Ref. [13] is extended to a longitudinal and a tilted 
vibration. 

The paper is organized as follows. In Sec. [Til the prob- 
lem is formulated: the governing equations and the di- 
mensionless parameters are introduced for the case of 
the longitudinal vibration. Also in this section, using 
the separation of the time scales, we split the nonlin- 
ear boundary value problem for the fluid flow into two 
coupled boundary value problems for the pulsatile and 
for the averaged flows. The pulsatile flow is analyzed 
in Sec. IIIII The averaged amplitude equation for the 
film height is obtained in Sec. IIVI using the solution of 
the pulsatile flow. The linear stability problem for the 
amplitude equation, the weakly nonlinear analysis and 
the numerical results on film dynamics are presented in 
Sec. |Vl Stability of the pulsatile fiow is demonstrated in 
Section IVII This stability translates into a stability of 
the averaged fiow and thus the averaged amplitude equa- 
tion is validated. In Sec. IVIII the analysis of the previous 
sections is generalized to the case of a tilted vibration. 
The conclusions are presented in Sec. IVIIIl 
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FIG. 1: Problem geometry: the longitudinal vibration. 



II. FORMULATION OF THE PROBLEM 

We consider a three-dimensional (3D) thin liquid film 
of the unperturbed height Hq on a planar, horizontal 
substrate. The Cartesian reference frame is chosen such 
that the x and y axes are in the substrate plane and the 
z axis is normal to the substrate (Fig.[T]). 

The substrate is subjected to the longitudinal har- 
monic vibration of the amplitude b and the frequency uj. 
We assume that the system is not confined in the horizon- 
tal directions or the vertical boundaries are motionless. 
Thus the substrate motion induces the fluid motion due 
to viscosity and sloshing modes are not excited. 

We assume the film height sufficiently small, so that 
the intermolecular interaction becomes important. In 
this paper, as in the preceding paper [l3l | we consider 
the van der Waals attractive potential. Generalization 
to other models of wetting interactions is straightfor- 
ward. Using same scalings as in Ref. [l^ [i.e., the units 
for the time, the length, the velocity and the pressure 
are Hq/v, Hq, v/Hq, p{v/HqY, respectively; u is the 
kinematic viscosity and p is the density of the liquid] we 
begin with the following dimensionless boundary-value 
problem: 



V • V 
V • Vv 



0, 



(la) 
(lb) 



V = BQ. sin ^ltex at z 
ht +v ■ V/i = w, 



0, 



{p + (p- CaK) n = n • T at z 



(2a) 

h{x,y,t). (2b) 



Here, v = (u, w) is the fluid velocity (where u is velocity 
in the substrate plane and w is the component normal to 
the substrate), p is the pressure in the liquid, T is the 
viscous stress tensor, h is the dimensionless height of the 
film, ex,z are the unit vectors directed along the x and 

z axes, respectively, n = {gz — V/i) / \J 1 -f {Wh f is the 
normal unit vector to the free surface, K = W ■ n is the 
mean curvature of the free surface, (j> — ~A/h^ [where 
A = A' / {dTrpiy^Ho) is the non-dimensional Hamaker con- 
stant], Ca = (7Hq/(pv^) is the capillary number (where 
a is the surface tension). Go — gaH^/v^ is the Galileo 
number, B = b/Ho is the non-dimensional amplitude, 
and fl = LoHg/iy is the non-dimensional frequency. 



We consider the nonlinear evolution of the large-scale 
perturbations. Proceeding exactly as in Ref. [l^, we first 
introduce a small parameter e, which is of the order of the 
ratio of the mean height Hq to the perturbation wave- 
length, i.e. e <§; 1 for long waves. Next, we introduce 
conventional stretched coordinates and the time 



X 



Y 



T 



assume large capillary number Ca — Ct~^ , and then sep- 
arate the pulsations depending on the "fast" time t = ^t 
and the averaged variables which depend on the "slow" 
time T. The detailed analysis of this procedure is pre- 
sented in it results in 



u = eU + U, w = e^W + (W, 
p — p + e^^p, h^h + eh, 



(3a) 
(3b) 



where all fields are 0(1) quantities. The pulsations (av- 
eraged variables) are marked by tildes (overbars). Sub- 
stitution of Eqs. ([3]) in Eqs. ^ and ([2]) gives two sets of 
equations and boundary conditions for the pulsational 
and the averaged parts of the velocity, pressure and 
height. 

(i) For the pulsations: 

= -V • U, ntSr = -Vp + tzz, (4a) 

Pz = 0, (4b) 

U = Bnainre^, = at Z = 0, (4c) 

nJir = -\J-Vh + W, 

\Jz = 0, p = at Z ^h, (4d) 

Hereafter V = {dx,dY,0) is a two-dimensional projec- 
tion of the gradient operator onto the X — Y plane, 
(ii) For the averaged parts: 



Wz ^ ~S7- U, Pz = -Go, 

Vzz = Vp+(U- VU + iyUz), 

U = W ^ at Z = 0, 

p = -(pzh)-<f>{h)-CV^h, 

JiT = -tS -Vh- (UV/i) + W - 

= -{IJzzh) atZ^h. 



{Wzh), 



(5a) 
(5b) 
(5c) 



(5d) 



In the set (O the angular brackets denote averaging with 
respect to the fast time r. The boundary conditions at 
the free surface have been shifted at the mean position 
h. Moreover, we neglect all terms of order e as they 
are unimportant for the further analysis. Note that the 
boundary value problem governing the oscillatory mo- 
tion, Eqs. (HI) is linear despite the finite intensity of the 
oscillatory motion, see Eq. ([5a|) . Also it can be seen that 
the set (jH) is decoupled from the set ([5]) and thus the 
solution of the former set can be immediately found. It 
is worth noting that B = 0(1) (in contrast to Ref. p^ . 
where the amplitude of the vibration has to be large in 
order to provide a finite intensity of the longitudinal mo- 
tion.) 
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FIG. 3: The value of = \l{Z = h)\ vs. Q. 

III. PULSATILE MOTION 



A. Analysis of the general case 

Here we assume stability of the pulsatile motion (see 
Sec. IVII for the proof) and determine the solution of 
Eqs. (g]). We seek the solution in the form 

U = Bne^Rc [I{X, Y, Z) exp {ir)] , (6a) 
W = BnRe[K{X,Y,Z)exp{iT)], (6b) 
p = 0, h = BRe[H{X,Y)exp{iT)]. (6c) 

Substitution of this ansatz in Eqs. (|4]) gives the set of 
equations and boundary conditions governing the ampli- 
tudes of the pulsations: 



Kz = -Ix, Izz + a'' I = 0, 
I = -i, K = a.t Z = 0, 
iH = K - Ihx, Iz ^0 a.t Z 



(7a) 
(7b) 
(7c) 



where = —iQ,. 

The solution of the boundary value problem ([7]) is: 



/ 
K 



ih 



.cos a{h — Z) 
cos ah 
1 — cos aZ 



X- 



cos^ ah 



H 



i-x 



cos^ ah 



(8a) 
(8b) 



Note that the amplitudes K, /, and H , generally speak- 
ing, depend on Y via h, but only the X-component of 
Vh is important for the pulsatile motion. 

Figure[5]shows Ir and li for various values oiQ,. (Hence 
and henceforth we use subscripts "r" and "i" to denote 
the real and the imaginary parts, respectively.) Figure [3] 
presents Im = \I {Z = H) \ at the mean position h. Plot- 
ting these figures we use the local frequency O = O/i^, 
which is determined through the local thickness of the 
layer. It is obvious that rapidly decays with increase 
of the vibration frequency. Note that I^hx is the ampli- 
tude of the surface deviation, \H\ (where, of course, hx 
is not apriori known). 



B. Limiting cases of oscillatory motion 

Case ^ 1. At large Q. the solution of the pulsatile 
motion, Eq. is small beyond the Stokes boundary 
layer adjacent to the substrate. Indeed, taking an ob- 
vious relation 2 cos a « exp(iQ;) = exp a* into account 
(the asterisk denotes the complex conjugation; here we 
assume that > 0), one immediately arrives at 

/ = -ze-"*^, (9) 
K = -2ie'''i^-'^'')hx, H = 4e~'^"''''hx. (10) 



As Eq. (|T0)) does not satisfy the no-slip condition, near 
the substrate this asymptotic formula has to be rewritten. 
Expanding 1 — cos aZ in a power series at small aZ in 
Eq. (|8b| . we obtain 



K = 2nZ^e-'^°''''hx at Z 



0. 



(11) 



Of course, Eqs. (fTU)) and (fTT|) do not match, since they 
are the opposite cases {VUz 1 and V^Z <ti 1, respec- 
tively) of the high-frequency approximation, ^ 1, for 
Eq. (fgbll . 

Since only the exponentially small terms are neglected 
in Eqs. ([^. pU]) these asymptotic expressions can be ex- 
tended even to moderate fl with high accuracy. For in- 
stance, the line 4 in Fig. [2] is indistinguishable from the 
curve corresponding to Eq. ([9^) a.t ft = 100 (|a| = 10). 

Case ri ^ 1. In the limit of small frequency the solu- 
tion of the pulsatile problem, Eqs. ^ reads: 



I = 

K = 
H 



-i-^Z {2h- Z) 



ZiZ^ ~AhZ^ 



in' 



8h-^ 



1 + !^ - I2h^ 



= 'IX 



hx[l- inJi^ - 



(12a) 
(12b) 

(12c) 
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Terms up to fi^ are held in the expansion of the gen- 
eral case solution, Eqs. This accuracy is needed to 
provide the averaged effects at low frequency. The ex- 
pression for /, Eq. (|12ap . explains the coincidence of the 
line 1 in Fig. [2l[b) with the line li = —1. Indeed, the dif- 
ference between these line is proportional to fl^ ^ 10^^, 
which cannot be seen on the scale of the figure. On the 
contrary, the real part [see Fig. [D^a)] is proportional to 
il and the corresponding variations are sufficient. 

Case fl arbitrary and h = const. In this limit it 
is clear from Eqs. (jSbp that the oscillatory flow is one- 
dimensional (ID) and there are no oscillations of the sur- 
face height. Thus in this case the flow is the oscillatory 
Couette flow generated by the vibration of the substrate 
in a layer with the free surface. Moreover, this flow differs 
from the well-known oscillatory Poiseuille flow: 



cos a{h — Z) 
cos ah 



(13) 



only in an additive constant. 

It is also important to note that Eqs. Q describe the 
conventional "Stokes layer", i.e. the ID flow forced by 
a high-frequency oscillation of the rigid plane in a semi- 
infinite space. 
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FIG. 4; Dependence of the coefficient Qi entering Eq. ([16 
on 7. 

linear equation for h: 



Ho 

Qi 
91 



V • ( -h'VUo 



-B^n^Qih^hx)^, (16a) 



'(l){h)~C\/^h + Goh. 
3(2gi -7.92) 



7^ (cosh 7 + 0087)" 
sinh 7 sin 7, 172 = sinh 7 cos 7 -|- sin 7 cosh 7 



(16b) 
(16c) 



IV. AVERAGED MOTION 
A. Analysis of the general case 

Using Eqs. the problem for averaged fields, 

Eqs. can be rewritten as follows (hereafter the over- 
bars are omitted): 



Pz 

Uzz 
U 

hx 

Uz 
P 



-Go, Wz^-y- U, 
1 



= Vp+ -B'Vt'RG{I*Ix +K*Iz)e.,, (14b) 
= = at Z = 0, 



(14a) 
(14b) 
(14c) 



-ii?2l]Rc(/*^i7)e„ 
~(b~CV^h at Z ^h. 



where 7 = \/2flh. This equation is the central equation of 
the paper. The first three terms at the right-hand side are 
the conventional terms resulting from van der Waals at- 
traction, capillarity and gravity. The term proportional 
new term resulting from the longitudinal 
vibration of the substrate. 

The dependence of Qi is given in Fig. [H Note that 
Qi > except for the narrow interval 5.268 < 7 < 8.507. 
There is also an infinite set of such intervals (the second 
one is at 11.69 < 7 < 14.85), but the corresponding 
absolute values of Qi are very small. 

It is clear that the longitudinal vibration along the X 
axis assigns the preferential direction {X) 'm X — Y plane. 
Thus the 0(2) group of symmetry, which is characteris- 
tic for the amplitude equations governing the thin films 
dynamics, is broken for Eqs. (|16p . However, the derived 
amplitude equation is still invariant under the transfor- 
mation X ~X. 



(14d) 



The evolutionary equation for h [the first equation in 
Eqs. (jl4d[l ] can be rewritten in the form 

hr^-W- / VdZ - -B^nRe {PH)^ a± Z ^ h. (15) 
Jo 2 

Analytical integration of this set of equations is per- 
formed in Appendix |^ It results in the following non- 



B. Limiting cases SI ^ 1 and Q. 

We first notice that the averaged effects vanish in the 
limiting case of high frequency, f2 3> 1. With expo- 
nentially small error, the oscillatory velocity given by 
Eqs. ([9]) and is ID and uniform along the X axis. 
Indeed, the only reason for the variations of the flow ve- 
locity stems from the variation of the film height /i, and 
thus this nonuniformity is exponentially small. As all 
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known mechanisms are based on the existence of a gradi- 
ent of the kinetic energy for the pulsations [51] , it is clear 
that the uniform oscillatory flux is unable to produce the 
averaged flow. Therefore, an averaged flow cannot be 
produced neither in the boundary layer nor in the core 
region. This conclusion agrees well with the asymptotics 
of Qi at large 7: 

6 

Qi ~ (sin 7 + cos 7) e ^. 

7 

Thus the only relevant limiting case is the case of low- 
frequency, O <C 1. Integration of the boundary- value 
problem for Eq. (jl4bp results in the following solution 
(see Appendix [Al for details): 

U = -lzi2h- Z)VUa - ^^U^^^hxe,, (17a) 

120 ^ ^ ' 

Then the following equation for the film height is ob- 
tained: 

/iT = V • Q/^'Vn^ - ^B^^^ {h^hx)j, , (18) 

which agrees with the expansion of Qi « 7^/30 in 
Eq. p6c)) at smaU Q,. 

It is seen in Eq. (|18p that at il ^ 1 the vibration 
impact is determined by the squared amplitude of the 
pulsatile acceleration. This conclusion is expected in view 
of the similar result for the vertical vibration. 



V. FILM DYNAMICS 

A. Linear stability analysis 

It is convenient to rescale the amplitude equation (jl6ap 
using X = y/C/'iAX, Y = ^C/2,AY, and T = 
CI{?,A^)f. [Recall that = -Ah-'^.] This leads to 



(19) 

where V = B^D? /2A and Go = Go/3A. 

Seeking the solution of Eq. in the form h = 1 + ^, 
where ^ is a small perturbation, one obtains: 



VQi{lo)^xx- (20) 



Here 70 = ^/h — \/2D,. 

For the normal perturbation ^ proportional to 

exp {ikxX -I- ikyY — AT^ the decay rate A is: 




k 1.6 0.8 
(c) 



K 1.6 



FIG. 5: Decay rates A for Go = 3.33 x 10"*, kv =0,kx= k. 
(a) n = 0.1. (b) Q^2,{c)fl = 20. (a) and (b): dashed, solid 
and dash-dotted lines correspond to V = 0, 200, 500, respec- 
tively; (c): dashed, solid and dash-dotted lines correspond to 
= 0, 2, 5, respectively. 



where k^ — kx+ky is the squared wavenumber. In terms 
of the original unsealed variables {kx, ky) 

A = i [Go - 3A + Gfc2] fc2 - lB^n^Q,{jo)k%. (22) 

Since Qi is positive except for the narrow intervals 
of 7 (see Fig. 3]), longitudinal vibration destabilizes the 
film beyond these intervals [see Fig. EJa) and Fig. Efb)]. 
Furthermore, one can readily see that the vibration does 
not impact the behavior of perturbations with kx = 
[longitudinal rolls). Thus stabilization of the film by ap- 
plication of the longitudinal vibration is possible only in 
2D systems, where there is no flow in the l"-direction, i.e. 
k = kx, ky — 0. Below we consider only the behavior 
of 2D perturbations, which are critical for the reasonable 
interval of frequencies. 

The typical pictures of the decay rate for this case are 
shown in Fig. [51 For stabilization [Fig. O^c)] one needs an 
extremely large frequency of the vibration. (Recall that 
r2 = 20 corresponds to n = 300 MHz for the water layer 
of the height lOOOl.) 

It is obvious that fco: 



fco = 1 + VQi (70) - Go 



(23) 



A 



(Go-l + P)P-yQi(7o)fcx, (21) 



solves an algebraic equation A = 0. For fc < fco an insta- 
bility takes place. 



B. Weakly- nonlinear analysis of 2D systems 

Since only the monotonic instability is present, as the 
analysis in Sec . I V Al confirms . we need to analyze branch- 
ing of stationary solutions. Based upon the results of 
Sec. IV Al we consider only the 2D system, i.e. h = h{X) 
for the stationary solution. Analyzing the periodic so- 
lutions of a given wavenumber k, we keep in mind that 
the obtained results are appropriate for the confined sys- 
tems of the length L = -k /k. Indeed, due to symmetry the 
boundary conditions hj^ — Q are imposed at A = 0, i for 
the periodic solution of given wavenumber k. The same 
boundary conditions should be set at the impermeable 
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boundaries X = 0, L for the confined system. In this case 
the spectrum of the wavenumbers for the perturbations 
is discrete and bounded from below: fc'"^ > k'^^^ = k, 
where A:*-"^ is the nth eigen wavenumbcr corresponding to 
(n + l)/2 wavelength confined into the horizontal length 
of the system. Thus for k > ko the longwave instability 
does not occur. Moreover, as it is shown below, even for 
k < kg the growth of perturbations does not necessarily 
lead to a rupture in a confined system. 

In the stationary case Eq. p9)) can be integrated once. 
Due to symmetry, the integration constant is set equal 
to zero. Thus we obtain 



2,/ll 



k^h 



V^~Go]h' = 0. (24) 



40 n 
V 

30 H 



20 



10 




1.5 



-r 

2 



— I 1 1 

2.5 ^ 3 



Here the primes denote the derivatives with respect to 
C = kX, i.e. the solution is assumed to be 27r-periodic in 

C- 

To study the weakly nonlinear behavior of the per- 
turbation, we expand the surface deflection h and the 
wavenumber k in powers of small 5: 



h^l + S^i+ S^^2 + k = ko + S^k2 



(25) 



Substituting these expansions in Eq. ((24|) we collect the 
terms of equal order in S. The first-order equation is 

LCi = kl^T +h + VQi (70) - Gol £,[ = 0. (26) 



Its solution has the form: 

~ a cos C, 



(27) 



whereas k^ is given by Eq. 
The second order equation is 



i6 = (4-M)aci, 



(28) 



where 



^1 ^ 7o 



2 d Qi (70) 



d-yo 70 



The solution of Eq. ^ is 



4 - VFi . 

6 = cos2C. 

^ 12fc2 ^ 



The third-order equation is: 



i6 = 

where 



(4 - fFi) ei6 - ^ (10 + VF2) 



2kokiS,'{ , 
(29) 



^ 7g Qi (70) 
^ 2 d7^ 70 



FIG. 6: The domain in the parameter space, where the su- 
percritical branching of the surface deviation takes place. 
Go = 0, but the curve is almost unchanged even for Go = 0.5. 



The solvability condition of this equation couples the cor- 
rection and the amplitude of the surface deviation as 
follows: 



k2 = 



10 + VF2 



a 

8ko' 



(30) 



If the term in the bracket is positive, then ^2 > and 
the subcritical bifurcation takes place. It is obvious that 
the branching solution is unstable in this case. Other- 
wise, the supercritical bifurcation occurs and the stable 
stationary solution corresponding to the deflected surface 
emerges. 

The curve separating these two regions in the plane 
{V, il) is shown in Fig. [S] It can be readily seen that the 
supercritical excitation exists only at the large enough 
values of V, when the destabilization effects are well pro- 
nounced. Nevertheless, this phenomenon is quite inter- 
esting and unexpected and requires an additional analy- 



C. Stationary periodic solutions 

To study stationary periodic solutions of a finite ampli- 
tude we integrate Eq. (|24p with the boundary conditions: 



/i' = at C = 0, 2tt, 



and the mass conservation condition: 

r2ir 

hdC = 27r. 





(31) 



(32) 



(Recall that for the equivalent confined system only half 
of the period should be taken, i.e. C < tt.) This boundary 
value problem was solved by the shooting method. The 
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FIG. 7: Variation of the maximal (the upper branches) and 
the minimal (the lower branches) thickness of the layer with k 
for V = 10 and = 1, 1.5, 2 - solid, dashed, and dotted lines, 
respectively. For the former case {Q = 1) the value of ko is 
marked. Stability of the corresponding solutions is indicated 
by arrows. 
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V \ 
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1 



1.8 2.2 2.6 3 k 3.4 , K k 



FIG. 8: Variation of the maximal (the upper branches) and 
the minimal (the lower branches) thickness of the layer with 
k. (a) - V = 30, f2 = 1, 1.5, 2 - solid, dashed, and dotted 
lines, respectively, (b) - sketch with marked fco and kc- solid 
lines correspond to stable solutions, dashed lines - to unstable 
ones. Domains of attraction are shown by arrows. 



numerical results are presented in Figs. [TlfTUl for Go = 
3.33 X 10-*. 

The amplitude curves h„i{k) are shown for ^ = 10 
and y = 30 in Figs. [7] and [51 respectively. These figures 
confirm the results of the weakly-nonlinear analysis: the 
supercritical bifurcation takes place for V — 30. 

In Fig. [7] the value of ko is marked for — 1 (solid 
line). For smaller k instability of the flat surface gives 
rise to a rupture. The stationary solution exists only for 
k > ko, i.e. all the curves in Fig. [7] represent unstable sub- 
critical solutions. The lower branches of the amplitude 
curves in Fig.[7]can be thought of as the boundaries of the 
domains of attraction for the equilibrium state {h = 1) 
in the framework of the evolutionary problem, Eq. (jl9[) 
(or the similar equation with C, = kX being introduced). 



FIG. 9: The shape of the free surface for V = 30, il = 2; 
fc = 3, 2.7, 2.56 - dotted, dashed, and solid lines, respectively. 
The solid line with smaller surface deformation corresponds to 
the stable state, another solid curve corresponds to unstable 
solutions as well as dotted and dashed ones. 

An initial distortion of the surface ho{C) with the trough 
deeper than hm necessarily leads to rupture, while an 
initial distortion with /io(C) > for any C, < C < 27r 
decays with time resulting in the equilibrium state at 
T ^ oo (see arrows in Fig. [7]). Of course, this interpreta- 
tion is not exact, as the whole variety of the initial states 
/io(C) is characterized by the only value, hm- Therefore, 
the domain of attraction has to be confined by a band 
of finite thickness. However, due to the fast growth of 
the attracting van der Waals potential with the decrease 
of h, hm is a perfect characteristic for the domain of at- 
traction and the band thickness is rather small. Our nu- 
merical tests based on the finite-difference computation 
of Eq. support this conclusion. 

For larger values of V the stable distorted surface is 
found within some interval, kc < k < ko, see the dashed 
and the dotted curves in Fig. ^a) and a schematic plot 
in Fig. [5]Jb). In the latter figure domains of attraction 
and stability properties of the obtained solutions are also 
demonstrated. Again, the lower branch hm{k) of the 
unstable solution can be thought of as the boundary of 
the domain of attraction: the rupture occurs for initial 
distortions with min/io(C) < hm', in the opposite case 
the initial perturbation decays and the stable branch is 
achieved (either the stable branch of solution for k < ko 
or = 1 for fc > fco). To the best of our knowledge, this 
problem provides the first example where a nonplanar 
free surface is stable in the presence of the attracting van 
der Waals potential. 

The interval of existence of a stable solution is quite 
small and for k — k^ a, saddle-node bifurcation takes 
place, see Fig. |Hl[b). Both stable (with smaller values 
of — 1|) and unstable (with larger values of — 1|) 
branches of the solution disappear at k = kc and do not 
exist at smaller k. 

The typical shapes of the surface along the dotted 
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FIG. 10: Variation of fc. with Q ior V ^ 20, 30, 40 - lines 
1-3, respectively. Here fc* = fco correspond to the pitchfork 
bifurcation (dashed lines) and k, — kc (solid lines) correspond 
to the saddle-node bifurcations. Dash-dotted lines show the 
locus of cusp points, where fco = fee; these lines correspond to 
ko{V) along the line shown in Fig. [S] 



curve in Fig. [H^a) is presented in Fig. [9l For the smaller 
value of the wavenumber, k — 2.56, bistability takes 
place, i.e. there are two stationary shapes, the stable 
and unstable one (with the smaller and larger surface 
deformation, respectively). 

The bifurcation lines ko{fl) and kd^) are shown in 
Fig. [ini The stable states with the deformed free surface 
exist between the corresponding solid and dashed lines. 
The direct pitchfork bifurcation takes place within this 
interval of il. Otherwise, the branching is subcritical. 

It should be noted that kc > kQ'^\ where kg'^^ = fco for 
V = 0, i.e. the critical wavenumber in absence of the 
vibration. This inequality ensures that rupture always 



occurs at fc < fcp'^'', i.e. there is no stable stationary states 
even with deformed free surface. Thus the longitudinal 
vibration cannot be applied for "nonlinear stabilization" , 
i.e. in order to produce stable distorted film with fc < 



VI. STABILITY OF THE TIME-PERIODIC 
MOTION (PULSATILE FLOW) 

A. Reduction of the stability problem 

The finite intensity of the time-periodic solution, 
Eqs. ([3]), (O, and ([5]), raises question of the solution sta- 
bility. In order to show stability or instability, we return 
to the unsealed Eqs. ([T]) and ^ and represent the veloc- 
ity, the pressure and the surface deflection in the form 

V = U + eWe^ + V, p = p + P, h^h + eh + T. (33) 

(Within this section we restore the bar for the averaged 
film height, h.) Here in accordance with Eqs. ^ the 
dominant parts of the unperturbed velocity are the pul- 
satile one, while for the pressure field and the layer height 
the averaged parts dominate over the oscillatory ones. 
Next, we linearize Eqs. (H]) and ^ with respect to the 



small perturbations V, P and T, assuming 2D perturba- 
tion V — {U,Q,W). (We will show that such perturba- 
tions are critical in the sense that they occur prior to the 
onset of 3D instability.) Thus we arrive at the following 
boundary value problem: 



Ut = 
Wt = 
u = 

p = 

Here we set 



(34a) 

-UoU^ ~ UUocc - WUo. - + V2;7,(34b) 
-UoW:,^UWo^-P,+V^W, (34c) 
W ^ at z = 0, (34d) 

0, Tt + UqT^ + Uh^^W at z^h. (34e) 



Uq = BilRe [I{z) exp (int)] 



(35) 



according to Eq. Besides, we have neglected all terms 
of order e, for example W and h. Thus with accuracy 
0(e) the base flow can be thought of as ID flow. 

Components of the viscous stress tensor are linearized 
as follows near the unperturbed surface, z = h(X,T): 



T, 



(0 



no • (Ti -I- 9,ToT) • no + 2no • To • ni (36) 

no-(Ti+a,ToT)-to 
+ni-To-to-fno-To-ti, (37) 



where Ti is calculated over U and W, Tq is the viscous 
stress tensor of the base flow [i.e. Tq^x = '^Uqx, Tqxz = 
Uoz, Tqzz = 0(e)], no and ni (to and ti) are the nor- 
mal (tangential) vector for the unperturbed free surface 
z — h and the correction due to its perturbation T, re- 
spectively: 



no 



ni 



hi 



1 - 



hi 



l + hl 



Finally, for the curvature we obtain: 



(l + hl 



3/2 



(1 + hi) 



5/2 ■ 



(38) 
(39) 
(40) 

(41) 



It is clear that the x-derivatives of the base state vari- 
ables are 0(e) quantities and thus they can be omitted 
from the leading-order analysis. Indeed, the base state 
depends on x only via the slow coordinate X, whereas 
the characteristic wavelength of perturbations is 0(1). 
For example, Uox = eUox- Also, for the same reason the 
derivatives with respect to the slow time T are omitted. 
This procedure is equivalent to the methods of frozen co- 
efficients: neglecting e-terms above we assume that the 
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base flow is ID and locally a;- independent. Moreover, the 
free surface is almost horizontal: according to Eqs. ([55]) 
one obtains no = and to = e^^ , which significantly sim- 
plifies the formulas. As a result, the stability problem for 
the pulsatile motion reads: 



Wt = -UoW^-P,+V^W, 
U = Ty = Oatz = 0, 



(42a) 
(42b) 
(42c) 
(42d) 



P 



c 



Tt = -f/oT^ + at z = /i 



T.. + 2(iy, -C/^T,)^ 



(42e) 



The primes denote the z-derivatives of Uq. Equa- 
tion (|14ap is taken into account to calculate entering 
Eq. dsn]). Recall that C = e^Ca. 

Large factor in normal stress balance [the first re- 
lation from Eqs. (|42ep ] results in T = 0, i.e. the surface 
is locally undeformable. This, in turn, means that the 
boundary value problem ([1^ is the stability problem for 
ID flow, periodic in time. This problem, the so-called 
time-dependent Orr-Sommerfeld problem, is analyzed in 
Sees. lyTBl and IVTC] 



B. Orr-Sommerfeld problem for the pulsatile flow 

It is well-known (31 that 2D perturbations are critical 
for either stationary or time-dependent Orr-Sommerfeld 
problem. This proves neglecting the y-component of V 
and allows to introduce in Eqs. (|42p a streamfunction ijj. 
Setting 

T = 0, C/ = -^., W^^^ 

and separating the x-coordinate by means of ip{x, z,t) = 
-0(2, t) exp (ikx) we arrive at 

D^i^t = -ik (UoD^iP - Uoijj^ + D^, (43a) 

^ = = at z = 0, (43b) 
^ = D'^ip = at z = /i, (43c) 

where D^ip = ip" — k'^if^. It is convenient to rescale the 
vertical coordinate and the time in such a way that h = 1. 
Keeping the same notations for the rescaled variables we 
obtain: 

D'^i^t = -ik (UoD^i; - C/^'?a) + D^, (44a) 

= ?/>' = at z = 0, (44b) 
■iP = 1)2^ = at z = 1, (44c) 



where 
Uo = 



RRe 



/o(z) exp I iftt 



.cos (5(1 — z) 



We again use the local oscillation frequency (cf. Sec. IIII|) 

ri = Qh^ , a — ah; 

the Reynolds number R is introduced as follows 

R = Bflh. (45) 

First, we show that any constant value Co can be added 
to the amplitude of the velocity oscillations Iq(z). It 
can be readily seen that Eqs. (|44|) are invariant under 
transformation 



t/j ip exp 



^Re (iCoe'^*) 



lo^Io + Co. (46) 



Because ip is transformed by the periodical in time factor, 
stability properties do not change under this transforma- 
tion. In particular, setting Co = « we reduce Iq to the 
velocity profile of the oscillatory Poiseuille flow, Eq. (E 



or 



lo 



cos q;(1 



cos a 



in terms of rescaled coordinate. 

Thus, the stability problem for the pulsatile flow is 
reduced to the Orr-Sommerfeld problem for the os- 
cillatory Poiseuille flow. In spite of the detailed studies 
of the modulated Poiseuille flow [H, [M 113, El, there 
have been no papers that deal with the particular case of 
the zero mean value of the pressure gradient. 

Intuitively, the flow in a finite layer is more stable than 
the flow in a semi-inflnite space, i.e. the Stokes layer, 
which is known to be stable [l3|- Nevertheless, we give 
some results on the stability of the flow in Sec. IVI C[ 
which conflrm the above guess. 



Stability problem for the oscillatory Poiseuille 
flow 



To solve the linear stability problem ([44]) for the time- 
dependent flow, we apply the following method. First, 
due to Floquet theorem the amplitude of streamfunction 
■0 can be represented in the form: 



-At 



*(z,i). 



(47) 



Here ^(z,t) is the periodic function of time with the 
period 2tt/(1, which thus can be expanded in a Fourier 
series as follows: 



+ 00 



*(z,t) ^*„(z)exp 



(48) 



cos a 



Substituting this ansatz into Eqs. (I44p . we arrive at the 
chain of coupled boundary value problems for the Fourier 
component, ^'„(z). Truncation of the series by replacing 
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FIG. 11: Floquet exponents A vs. k for Re = 50000. Three 
lower branches of the spectrum (a)-(c) for Cl = 0.2, 1, 5 - 
dotted, sohd, and dashed hnes, respectively. 



the upper and lower limits of summation with N{—N), 
respectively, leads to the set of 4(2A^+ 1) ordinary differ- 
ential equations. This set has been solved by the shoot- 
ing method. We followed several lower branches of the 
spectrum; well-pronounced stabilization of the flow with 
increase of R was found for all the branches. 

The example of the computations is presented in 
Fig. [TTl where the three lower branches of the Floquet 
exponent A are shown. The curves are obtained with 
N — 10; the changes caused by larger N cannot be 
seen on the scale of the figure. Therefore, the oscilla- 
tory Poiseuille flow is shown to be stable even for finite 
values of fl. 



VII. VIBRATION ALONG THE 
NEAR- VERTICAL AXIS 

General notes and the analysis of the pulsatile 
motion 



In this section we brie fly g eneralize the previous analy- 
sis (Sees. HiniVl and Ref. ^^3]) to the case of the vibration 
along the axis that is tilted at a certain angle (3 to the 
substrate (see Fig. [T2)) . It is obvious that at the finite 
values of /3 the normal component of the acceleration is 
unimportant. Indeed, it is shown in Ref. [isj that the 
impact of vertical vibration is finite only at large am- 
plitudes: B = 0(e~^), while the horizontal vibration 
becomes essential even at finite B. Therefore, the longi- 
tudinal component of the vibration velocity is determi- 
native. 



\ 


4. 








H{x,y,t) 



FIG. 12: Problem geometry: the tilted vibration. 



The only case where the effect of the normal acceler- 
ation can compete with the one due to longitudinal mo- 
tion is the large-amplitude almost vertical vibration, i.e. 
/3 ^ 0. This case is important for applications because 
in experiments it is difficult to ensure the absolutely ver- 
tical axis of the oscillatory motion - the horizontal com- 
ponents of accelerations occur inevitably. 

Thus we assume (3 small, i.e. cos /3 « 1 and the near 
vertical motion of the substrate according to the law z — 
zq — B cos T (in the laboratory reference frame) . The 
fluid motion in the reference frame moving vertically with 
the substrate is governed by the Eqs. (H]) and the 
replacements are (i) the gravity modulation: 



G(<) = Ga + Bfl^cosT 



(49) 



and (ii) the longitudinal velocity is now (3Bfl sin flte^ in 
Eq. dial). 

Assuming large vibration frequency B we introduce the 
rescaled amplitude of the vibration b = e^^B and the 
rescaled angle of the vibration /i = e~^/3. Thus the am- 
plitude of the longitudinal motion of the substrate is fib 
and the limiting case fj, = corresponds to the transversal 
vibration, i.e. the results of Ref. [1J| are reproduced. For 
the longitudinal vibration it is necessary to set ^ 1, 
6^1, while keeping the product, B = finite in order 
to obtain the corresponding formulas from Sees. IllUfVl 

Representing all fields as the sums of the pulsatile and 
averaged parts according to Eqs. ([3]) we arrive at the 
following equations for the pulsations: 

Wz = -S7- U, niJr - -Vp + iJzz, (50a) 

pz = -bn'^cosT, (50b) 

U = bmnsmr, T4^ = at Z = 0, (50c) 

nJir = -tS-Vh + W, 

tSz = 0, p = &t Z = h. (50d) 

Here m = iie^- The averaged motion is described by 
Eqs. p. 

Representing the solution of the boundary value prob- 
lem ([50)) in the form 

p = 5f7Re [g(X, F, Z) exp (zt)] , (51a) 

U = 6f7Re [I(A:, r, Z) exp (ir)] , (51b) 

W = bnRe[K{X,Y,Z)exp{iT)], (51c) 

h = 5Re [i7(A:, r) exp (ir)] . (51d) 
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[cf. Eqs. (20) in [T^ and Eqs. ([6])] and solving the equa- 
tions for the amphtudes g, I, if, and H wc obtain (here- 
after the bar over h is omitted again): 



q 

I 

K 



n{h- z), 



iVh ~ iM 



+iVh ■ M 



cos a{h — Z) 

cos ah ' 

sin a (h ^ Z) — sin ah 

+ — ^ — 4 

a cos ah 
1 — cos aZ 



cos^ ah 



H = -V • [/ {ah) hS/h] 
where M = V/i -t- m and 



m • V/i 

cos^ ah ' 



(52a) 
(52b) 

(52c) 
(52d) 

(53) 



Due to the hnearity of Eqs. ([50)1 the oscillatory flow is a 
superposition of the motions generated by a vertical [see 
Eqs. (42) in [ll]] and a horizontal [Eqs. ^] vibration. 

Referring to Sec. IIII Bi as well as to Sec. Ill of Ref. |l3| . 
we do not present the limiting cases of the solution ([5^ . 
Milestones of the stability analysis for the pulsatile ve- 
locity are given in Appendix [Bl 

Generally speaking the solution (|52p remains valid even 
for the complex-valued /i and m. This permits to con- 
sider the vertical and horizontal vibration to be out-of- 
phase, but leads to the cumbersome equation for the av- 
eraged fields. Thus, hereafter we assume real m. 



B. Amplitude equation and limiting cases 

Solving the averaged problem using the procedure sim- 
ilar to the one in Appendix[Xl we arrive at the amplitude 
equation for the averaged height: 



hT 

n 
Q 



Ho 
h' 



i.3vn) - 



1. 



b^n^v ■ Q, 



2 

Q2iV^hM 



1 



(54a) 
(54b) 



g22M + -m) • VV/i 



+h^QiMM-\/h, 
RcH = -frhV'^h - {Vhf + Q3M • V/i. 



(54c) 
(54d) 



Here Hq and Qi are given by Eqs. (jl6bp and (|16cp . re- 
spectively, 



?2i = 6gi-2g2, Q22 = 5gi - 92 
cosh 7 cos 7 -|- 1 



(cosh 7 -I- cos 7) 
fr = Ref (ah) = 1 



sinh 7 + sin 7 



7 (cosh 7 4- cos 7) 



sinh 7 — sin 7 
7^ (cosh 7 + cos 7) 



cosh 7 — cos 7 



7^ (cosh 7 + cos 7) 



(55) 
(56) 

(57) 




FIG. 13: Left axis - the dependence of the coefficients Q22 
(solid line) and Q3 (dashed line) in Eq. (|54|) . on 7. Right axis 
- Q21 (7). 



The variation of the coefficients (32i, Q22, and Q3 with 7 
is shown in Fig. I13i whereas fr can be found in Fig. 6(a) 
of Ref. 13]. 

For 2D case, when h depends on X and T only, the 
amplitude equation reduces to 



^h'Uxj^-^b'n^Qx, 



(59a) 



Q ^ Qih)h^{hx + f.ifh 



X 



+Q2 (7) h^hxx {hx +f^) + ^f^h^hxx (59b) 

Here Q2 = Q21 + Q22, cf. Eq. (31b) in Ref. [H. 

Equations (|54p coincide with Eq. at \fi\ ^ 1, jib = 
B (the longitudinal vibration). Equations ((59l) for /i ~ 
(the vertical vibration) coincide with Eqs. (30) in 
Ref. [3], but the corresponding 3D analogue reads 



hr 

Q 



V • ( ^/i^VH 



Qih^ [VhyVh 



(60a) 



'h^ [Q2iV2/iV/i + Q22V/1 • VV/i] , (60b) 



M must be replaced with Vh in the definition of H, 
Eq. (|54dp . It is clear that Q differs from the one defined 
by Eq. (43c) in Ref. [13]. This contradiction is caused by 
the calculation mistake in Ref. p^ . 

It is clearly seen that the amplitude equation (|54p con- 
tains several cross terms, which are linear with respect to 
fi and proportional to the first derivative of h with respect 
to X. These terms remove the degeneracy with respect 
to the replacement X — > —X. Only the invariance under 
the simultaneous transformation X —X and fj, —^ — /i 
holds. Thus the presence of the vertical vibration makes 
different the motion along the positive and negative di- 
rections of the X axis. Indeed, the X-component of the 
oscillatory velocity is in-phase with the vertical compo- 
nent, whereas the projection of the pulsatile velocity on 
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the —X axis is in counter-phase. This phase shift results 
in the difference after averaging. 

The corresponding Hmits for the general case of the 
tilted vibration give the following amplitude equations: 

(i) low frequency (fi ^ 1): 



inary parts of the decay rate we arrive at 



Q; 

n 



iv • (/i^Vn) + h^n'^V ■ Qi, (61a) 

o 

315 



y^hM + 9M-\7\7h) 



—MM ■ V/i, 
15 

no - . {h^M) 

15 



(61b) 
(61c) 



For the vertical vibration this leads to Eq. (|61ap with 



Qi 



n 



^ (V^hVh + 9Vh-VVh) 
315 

-—{VhfVh, 



15 



Tin 



(62a) 
(62b) 



instead of Eq. (44) in [13]. For the horizontal vibration 
Eq. (fT5|) is reproduced. 

(ii) high frequency {il ^ 1): 



n 



3 
Ho 



V • (hVh) . 



(63a) 
(63b) 



It is obvious that the longitudinal component of the vi- 
bration has no impact in this limiting case (see the expla- 
nation in Sec. IIVB|) . Therefore, the obtained expression 
coincides with the corresponding equations obtained in 
Ref. [HI as well as Eqs. (46) and (47) in [11]. 



1 



Go-3A + k^{C 



-fr 



2 

b^fl^ 



fi Qik 



sinh 7 — sin 7 
3 (cosh 7 -I- cos 7) 



Q7 



(65) 
(66) 



It can be readily seen that the impact of tilted vibra- 
tion on the real part of the decay rate, Ar, is the su- 
perposition of impacts from the vertical and horizontal 
vibration [cf. Eq. (72) in [Ts] for the former case and 
Eq. ((22| for the latter one]. The additional terms, linear 
with respect to /i, provide only the imaginary part of A. 

Thus the vertical component of the vibration is unim- 
portant for the longwave perturbations (with /c — > 0), and 
in this case the tilted vibration, similar to the horizontal 
one, decreases the stability threshold (unless Qi < 0, see 
the discussion of Eq. (|^ ). 

Again, the 2D perturbations [ky = 0, kx = k) are 
critical for Qi > and longitudinal rolls {kx = 0) are 
critical for narrow intervals of f2, where Qi < 0. However, 
the latter case seems unrealistic as very high frequencies 
are needed. 

In confined systems with a discrete spectrum of k the 
competition of the stabilizing effect of the vertical vi- 
bration and the destabilizing effect of the horizontal one 
takes place. 

It should be noted also that the emergence of the imag- 
inary part of A is the indicator of the averaged transport 
in the system. Indeed, the perturbations are stationary 
in a reference frame moving along the X axis with a con- 
stant velocity Xijkx- However, such a longitudinal drag 
is not limited by the transport of perturbations: any ad- 
mixture can be spread over the system by means of the 
tilted vibration. Thus the tilted vibration seems to be 
the novel way to transport microparticles or molecules, 
which is important in many microfluidic applications. 



C. Stability analysis of the flat surface 

Representing h in the form 1 + C and linearizing 
Eqs. (|54p with respect to a small perturbation ^ results 
in 



1. 



(Go-3A)e- 



C 



b'il' f3Q2 + l~Qs 



f^y^^x + m'QiCxx j64) 



Here 7 must be replaced by 70 in Qj, j = 1,2,3 and 
fr, because all the coefficients are calculated for the un- 
perturbed state, Hq = 1. Substituting ^ proportional to 
exp {ikxX + ikyY — AT) and separating real and imag- 



VIII. SUMMARY 

We consider a thin liquid film on a planar horizontal 
substrate subjected to a high frequency vibration. In the 
absence of a vibration, the van der Waals attraction to 
the substrate destabilizes the film and causes its dewet- 
ting. In contrast to conventional averaging method, we 
assume that the period of the vibration is comparable to 
the time of viscous relaxation of perturbations across the 
layer. This allows us to apply the averaging method to 
the ultra-thin films. Such analysis was first developed in 
Ref. [ll], where the vertical vibration is considered and 
is shown to enhance film stability. 

This work is a natural extension of Ref. |13|] . We con- 
sider, separately, the longitudinal and the tilted vibra- 
tion. In the former case the finite amplitude of the vibra- 
tion results in destabilization of the layer. There is also 
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a sequence of narrow intervals of the vibration frequency, 
where stabihzation occurs in the two-dimensional prob- 
lem. However, the frequency must be very high (at least 
300MHz for a water layer of the thickness 1000 A). 

Another effect of the longitudinal vibration is the emer- 
gence of the supercritical branching at the sufhciently 
high intensity of the vibration. In this case the deformed 
free surface becomes stable, i.e. the instability of the flat 
surface does not necessarily lead to a rupture. 

For the tilted vibration the longitudinal (destabiliz- 
ing) component of the pulsatile velocity is shown to be 
dominant. The only case, where the competition of the 
vertical and the horizontal vibration occurs, is the almost 
vertical vibration of large amplitude. For this case the 
averaging procedure is carried out and the corresponding 
amplitude equation is obtained. This analysis allows to 
correct the three-dimensional generalization of the am- 
plitude equation for the vertical vibration [Tsj. 

Linear stability analysis in the framework of the am- 
plitude equation indicates that both destabilization and 
stabilization of the flat surface are possible in the case 
of tilted vibration. Stabilization takes place only in the 
confined systems, when the spectrum of perturbations is 
discrete and bounded from below. 

Besides, the small perturbations are oscillatory, i.e. the 
drag takes place for the tilted vibration. This property 
can be very important for many microfluidic applications 
since some admixtures can be transported in the same 
manner as the perturbations. 
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APPENDIX A: SOLUTION OF THE AVERAGED 
PROBLEM. LONGITUDINAL VIBRATION 

The solution of the averaged boundary value problem, 
Eqs. (HH), is sought in the form: 

p = -GoZ + Ho, (Al) 
U = U, + ^i7„e„ (A2) 
where Hq is given by Eq. (|16b[) . 

Uc = iz(z-2/i)vn 

is the "conventional" part of the longitudinal velocity, 
and its "vibrational" part is determined by the boundary 
value problem 

C/^' ^ Reiir^+Ki;), (A3a) 
Uv = Oat Z = 0, (A3b) 

K = -^^<={Hrzz) atZ^h. (A3c) 



It can be readily shown that 

Uy = Rc / KI*dZ - / d£,dx / |/pd77 (A4) 
Jo Jo J^ 

is the solution of Eqs. (|A3p . 

Substituting Eq. (|A2[) with both known parts, Uc and 
Uv, into Eq. (fT5|) one can obtain Eq. ((T6)) after straight- 
forward (but cumbersome) integration. 

APPENDIX B: STABILITY OF THE PULSATILE 
FLOW. TILTED VIBRATION 

Here we only briefly analyze the stability of the pul- 
satile flow ([52)1 as this analysis is quite similar to one 
given in Sec. VII of Ref. 13] and Sec. IVII of the present 
paper. 

In order to study the behavior of short- wavelength per- 
turbations, we again return to the unsealed system of 
equations ([T]) and ([2]) with the above-mentioned changes: 
(i) and (ii), see Sec. I VII Al The perturbed fields are 

V = U + eWe^ + V, p = e-^p + p + P, h = h + eh + T. 

(Bl) 

The only difference of Eqs. (jBl|) from Eqs. ([33]) is the 
presence of the oscillatory part of the pressure, according 
to Eq. (|52ap . In Eq. (|Bip the oscillatory parts are given 
by Eqs. ((52)) . Below we again omit the bar over h. 

Repeating the same procedure of the method of frozen 
coefficients, we arrive at the following set of equations 
governing the small-amplitude perturbations: 



w = 


-V-U, 


(B2a) 


Ut = 


-Uo • VU - 


WV'a - VP + V^U + U"(B2b) 


Wt - 


-Uo • VW - 


- P' + \7'^W + W", (B2c) 


u = 


W = at z - 


= 0, (B2d) 




-Uo • VT + 


W, U' + VW = -Ujj'T, 


P = 


fbn^ ^ 

1 cos ilt 


+ g + ^)t-5v^t 



-1-2 {W, - U;, • VT) at z^h. (B2e) 



Here the primes denote the z-derivatives and the base 
oscillatory velocity is 

Uo = bflRe [I(z) exp (int)] . (B3) 

This boundary value problem is a straightforward 3D ex- 
tension of Eqs. (|42|) . with two exceptions. Firstly, there 
is an additional term at the right-hand side of Eq. (|B2e[) . 
which comes from —p^T. Secondly, the base oscillatory 
velocity is 2D. 

As it is shown in Ref. [13] [see Eqs. (52) there] , the pres- 
ence of the two asymptotically large terms in Eq. (|B2e[) 
allows to split the stability problem into two different 
problems: (i) the Faraday instability and (ii) the insta- 
bility of the oscillatory flow with the undeformable free 
surface. 
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(i) The solution is based on the results by Mancebo 
and Vega Q. This analysis gives Eq. (53) of Ref. 13]: 



MV 



(B4) 



where ^mv is the function given in Fig. 4 of Ref. [8|. 

(ii) It is easy to see that the stability problem for 2D 
base flow (jB3[) can be reduced to Eqs. (jH]). Indeed, after 
the separation of the longitudinal coordinates according 
to f{x, y, z, t) = f{z, t) exp(ik-R) and the obvious trans- 
formation 

f ^ f exp [—ib cos fiik • Vh] , 
the base velocity (|B3|) becomes locally ID: 



Uo 



-fer^MRe 



xosajh- Z) 
cos an 



(B5) 



Here / is each of the fields {P,W,XJ - 
and R are 2D vectors in plane x — y. 



(t/,y),T}, andk 



Thus, we again exclude the constant part of the ve- 
locity, i.e. the first term in I by the periodical in time 
transformation similar to Eq. ((46|l (see Sec. I VI C|l . Now 
one can redirect the local x axis along the vector M and 
introduce the streamfunction in order to obtain exactly 
the same Orr-Sommerfeld problem as Eq. (|43)) . It has 
been shown in Sec. IVI Cl that the latter flow is stable. 
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